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INTRODUCTION 

Since  its  discovery  by  Powell  and  Stetson  in  1965,  time  average  holo- 
graphic interferometry  has  developed  into  an  extremely  useful  tool  for 
studying  the  resonant  mode  shapes  and  frequencies  of  vibrating  bodies. 

The  method  has  proven  to  be  particularly  useful  in  the  turbine  engine 
industry  where  the  vibration  properties  of  turbomachinery  blading  must  be 
known  in  order  to  insure  that  dangerous  blade  resonances  do  not  coincide 
with  operating  engine  speeds.  While  the  mode  shape  and  natural  frequency 
can  be  directly  determined  by  time  average  holographic  interferometry, 
knowledge  of  the  corresponding  strain  distribution  for  each  vibration  mode 
shape  is  also  desirable.  Computation  of  the  surface  strain  has  presented 
some  problems  in  the  past  because  of  the  error  magnification  that  results 
in  using  inexact  experimental  data  to  compute  the  second  derivative. 
Various  approaches  have  been  proposed  in  the  literature  for  attacking  this 
problem.  Waters,  Aas,  and  Erf  (Ref.  1)  determined  the  bending  strain  of 
the  first  bending  vibration  mode  of  a turbine  blade  using  a "smooth  curve" 
approach  in  order  to  smooth  out  errors  caused  by  differentiation.  Taylor 
and  Brandt  (Ref.  2)  carried  out  a very  good  error  analysis  of  strain  com- 
putations based  on  cubic  spline  functions.  Recently,  Dandliker,  Ineichen, 
and  Mastner  (Ref.  3)  reported  on  an  experimental  technique  for  measuring 
the  interference  phase  to  within  .3°  at  any  point  on  the  displacement 
fringe  pattern;  thus,  improving  the  accuracy  of  the  holographic  displace- 
ment data  by  two  orders  of  magnitude. 

In  the  vein  of  the  first  two  references  cited,  the  present  study  uses 
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standard  time  average  holographic  fringes  as  the  raw  data  for  determining 
the  bending  strain  of  a plate-like  structure,  e.g.,  turbine  engine 
blading.  The  holographic  fringes  are  used  as  the  raw  data  to  generate 
a mathematically  continuous  series  approximation  of  a plate-like  struc- 
ture’s normal  displacement  where  the  terms  of  the  series  approximation  are 
clamped-free  or  free-free  eigenfunctions  of  a simple  beam.  Since  the 
series  is  continuous,  it  can  be  differentiated  twice  to  obtain  the 
curvature  or  bending  strain  resulting  from  one  of  the  plate's  vibration 
mode  shapes.  The  use  of  beam  functions  in  the  series  approximation  has 
the  advantage  that  the  plate's  geometric  boundary  conditions  are  identi- 
cally satisfied  and  its  natural  boundary  conditions,  satisfying  equili- 
brium, are  approximately  satisfied. 

What  follows  is  a development  of  the  theory  upon  which  the  method 
is  based,  a description  of  the  resulting  interactive  computer  program, 

.and  some  examples  of  its  application  to  some  plate  vibration  modes.  The 
results  are  compared  with  bending  strain  values  obtained  using  finite 
element  analysis. 
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SECTION  II 
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SUMMARY  OF  RESULTS 


A technique  for  computing  the  bending  strain  resulting  from  the 
resonant  modal  deformation  of  vibrating  plate-like  structures  is  de- 
scribed. Interferometric  fringes  obtained  by  time  average  holography 
are  used  as  the  basis  for  generating  a mathematically  continuous  series 
approximation  of  a plate-like  structure's  normal  displacement.  The 
terms  of  the  series  consist  of  the  clamped-free  or  free-free  eigen- 
functions of  a simple  beam.  The  bending  strain  is  then  obtained  by 
computing  the  second  derivative  of  the  displacement  series. 

The  coefficients  of  the  displacement  series  terms  are  computed  for 
a given  segment  (length  or  width)  of  a cantilevered  plate-like  structure 
based  upon  the  holographic  fringe  values  lying  along  the  same  plate 
segment.  A linear  least  squares  solution  routine  is  used  to  solve 
for  p series  coefficients,  called  modal  weighting  coefficients,  in  terms 
of  k normal  displacement  values  obtained  from  k holographic  fringe  values 
where  k^^p.  A "best  fit"  solution  is  thus  obtained  for  the  plate  dis- 
placement. This  least  squares  approach  in  conjunction  with  the  fact 
the  beam  series  functions  exactly  satisfy  the  plate's  geometric  boundary 
conditions  and  approximately  satisfy  the  plate's  natural  boundary  condi- 
tions, results  in  a displacement  series  that  yields  quite  accurate  dis- 
placement and  bending  strain  values. 

The  technique  described  above  has  been  programmed  for  use  on  a 
Tektronix  4010  interactive  computer  terminal.  To  use  the  program,  called 
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HOLOCURVE,  the  user  specifies  the  particular  chordwise  (free-free 
boundary  conditions)  or  spanwise  (clamped-f ree  boundary  conditions) 
plate  segment  to  be  studied,  the  number  of  terms  to  be  used  in  the  series 
approximation,  and  the  holographic  fringe  values  along  the  plate  segment. 

A listing  and  plot  of  the  plate's  normal  displacement  and  bending  strain 
is  generated  by  HOLOCURVE. 

The  accuracy  and  effectiveness  of  the  technique  used  by  HOLOCURVE 
is  checked  by  determining  the  displacement  and  bending  strain  at  selected 
segment  locations  for  three  different  modes  of  vibration  of  a cantilevered 
plate.  The  results  are  compared  to  those  of  an  eigenvalue  analysis  carried 
out  on  a finite  element  model  of  the  plate  using  the  finite  element  com- 
puter program,  NASTRAN.  The  HOLOCURVE  results  are  in  excellent  agreement 
with  the  finite  element  analysis  except  for  cases  where  the  bending  strain 
is  essentially  zero  as  in  the  case  of  chordwise  segment  for  a torsional 
mode  shape.  In  cases  such  as  this,  HOLOCURVE  Issues  a warning  message  to 
the  user. 
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SECTION  III 


THEORETICAL  DEVELOPMENT 

3 . 1 Series  Approximation  of  the  Plate  Deformation' 

In  the  spirit  of  Rayleigh's  method  (Ref.  4),  Introduced  one  hundred 
years  ago  in  his  treatise,  "The  Theory  of  Sound",  one  can  approximate  the 
normal  deformation  of  a freely  vibrating  rectangular  plate  by  the  natural 
mode  shapes  of  beams  whose  boundary  conditions  are  similar  to  those  of  a 
plate  along  its  respective  edges.  Thus,  for  a cantilevered  plate 
(Figure  3.1),  having  one  edge  fixed  and  three  edges  free,  one  can  model  the 
plate's  modal  deformation  in  the  clamped-free  direction  (parallel  to  the 
X axis  in  Figure  3.1)  by  a similar  clamped-free  beam  mode  shape.  Similarly, 
the  plate's  modal  deformation  in  its  free-free  direction  (parallel  to  the 
y axis  in  Figure  3.1)  can  be  modeled  by  suitably  chosen  natural  mode  shapes 
of  a free-free  beam.  The  boundary  conditions  of  the  plate  and  the  corres- 
ponding clamped-free  and  free-free  beams  are  similar  in  the  sense  that  the 
geometric  boundary  conditions  (deflections,  rotation,  and  slope)  are  exactly 
satisfied  while  the  natural  boundary  conditions  (equilibrium  conditions)  are 
only  approximately  satisfied  (Ref.  5).  The  ramifications  of  only  approxi- 
mately satisfying  the  natural  boundary  conditiors  will  be  discussed  later. 

Referring  to  Figure  3.1,  let  the  plate's  normal  deflection  with  respect 
to  the  spanwlse  direction,  x,  (clamped-free  boundary  conditions)  be  repre- 
sented by 

P 

Wjj(x)  = AiXj:(x)  (3,1) 

where  i - mode  shape  number 

Xi  - 1-th  natural  mode  shape  for  a clamped-free  beam 
Ai  - modal  weighting  coefficient  for  mode  i. 


Figure  3.1  - Coordinate  System  for  Cantilever  Plate 


i>6 


Referring  once  more  to  Figure  3.1,  let  the  plate's  normal  deflec- 
tion with  respect  to  the  chordwise  direction,  y,  (free-free  boundary 
conditions)  be  represented  by: 

q 

Wy(y)  = BiYi^(y)  (3 

where  i - mode  shape  number 

- i-th  natural  mode  shape  for  a free-free  beam 
Bi  - modal  weighting  coefficient  for  mode  i 


For  a beam  of  length,  b,  having  free  boundary  conditions  at  y=0  and  y*b, 
the  natural  mode  shapes  are  given  by  the  following  expressions  (Ref.  6): 


Y^  = 1 (3.4) 

Y^  = (1  - -^)  (3.5) 


Yj^(y)  = coshli^  + cos.1^  - Sj^fsinh-^^^  4. 

b b b b J 

for  i = 3,  4,  5 ... 


(3.6) 


where  - eigenvalue  for  the  i-th  mode  shape  of  a free-free  beam 

- constant  dependent  upon  the  mode  number  and  boundary 
conditions. 


The  eigenvalue,  is  related  to  the  free-free  beam  properties 
in  the  same  way  as  was  for  a clamped-free  beam.  Table  3.2  gives 
the  eigenvalues,  Yi»  and  constant,  6i,  for  a free-free  beam. 
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TABLE  3.2 


i 

Free-Free  Beam  Constants  (Ref.  6) 

^i 

^i 

3 

4.7300 

.9825 

4 

7.8532 

1.0007 

5 

10.9956 

.9999 

6 

14.1371 

1.0000 

N^7 

(2N-3)it 

1.0 

2 

Note  that  the  expressions  given  by  equations  (3.4)  and  (3.5)  represent 
rigid  body  modes  of  the  free-free  beam.  Equation  (3.4)  is  a normalized  trans- 
lation of  the  beam  in  the  z direction  and  equation  (3.5)  represents  a rota- 
tion of  the  free-free  beam  about  its  midpoint,  the  vector  of  rotation  being 
parallel  to  the  x axis. 

The  modal  functions  given  by  equation  (3.2)  for  the  clamped-free  beam 
and  by  equations  (3.4)  - (3.6)  for  the  free-free  beam  are  continuous  with 
respect  to  each  function's  spatial  variable  and,  hence,  can  be  differentia- 
ted twice.  Thus,  from  equation  (3.1),  the  plate  curvature  in  the  x dlrec- 


similarly,  the  plate  curvature  In  the  y direction  can  be  repre- 
sented by 


Wy  ^ ®i^i 


1 = 1,  2,  3,  ...  q 


(3.9) 


y 1=1 

where  the  second  derivative  of  the  1-th  free-free  mode  Is  given  by 

for  1 = 1,  2 (3.10) 

Yl\2| 


(y)  = 0, 


(y)  = (tt) 


y Y 

— - cos 


— - 3 sinh - sin  ^ )1  (3.11) 

b b b ‘'j 

for  1 = 3,  4,  5,  . . . 


Neglecting  in-plane  stretching,  equations  (3.7)  and  (3.9)  can  be 
used  to  determine  the  bending  strains  in  the  plate  caused  by  the  out  of 
plane  deformation.  For  a plate  undergoing  pure  bending,  the  strains 
can  be  obtained  from  (Ref.  5): 


e = -zW  (x) 

X X 


Gy  = -zWy  (y) 


(3.12) 

(3.13) 


where  z is  the  perpendicular  distance  from  the  plate's  neutral  surface 
to  the  point  of  interest  on  the  plate's  surface.  Note  that  for  a 
twisted  plate-like  structure  such  as  a compressor  blade,  the  distance, 
z,  does  not  necessarily  equal  half  the  plate  thickness.  Using  equations 
(3.7)  and  (3.9),  equations  (3.12)  and  (3.13)  can  be  written  as: 

P 


e = -z  Z AjiX^  (x) 
X 1=1 


e„  = -z  I BiY,  (y) 
^ i:-l 


(3.14) 


(3.15) 


1 
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Looking  at  equations  (3,14)  and  (3.15),  we  see  that  the  normal 
strain  with  respect  to  x (y  held  constant,  i.e.,  constant  chord 
location  for  a blade)  or  the  normal  strain  with  respect  to  y (x  held 
constant,  i.e.,  constant  span  location  for  a blade)  can  be  computed 
provided  that  one  can  arrive  at  suitable  values  for  the  A^’s  and  B^'s. 
What  follows  is  a way  of  accurately  computing  the  Ai's  and  Bi's  for  a 
specified  natural  mode  of  vibration  of  a plate-like  structure  by  utili- 
zing the  interference  fringes  from  time  average  laser  holography. 


The  weighting  coefficients,  Ai  and  Bi,  that  appear  in  equations  (3.1) 
and  (3.3)  for  approximating  the  plate  deformation  and  also  in  equations 
(3.7)  and  (3.9)  for  approximating  the  plate  curvature  can  be  determined 
using  a least  squares  technique  based  on  discrete  displacement  data  points 
obtained  by  laser  holographic  Interferometry.  For  a plate-like  structure 
vibrating  in  one  of  its  natural  modes  of  vibration,  time  average  holographic 


interferometry  can  be  used  to  determine  the  modal  deformation  of  the  struc- 


ture (Ref,  7),  The  modal  deformation  is  characterized  by  a "contour  map" 
of  fringes  appearing  on  the  surface  of  the  plate-like  structure  when  the 
time  average  hologram  is  reconstructed  with  laser  light.  A photograph  of 
one  such  reconstruction  is  shown  in  Figure  3.2  which  shows  a turbine 
blade  vibrating  in  one  of  its  natural  modes  of  vibration.  The  fringes 
on  the  blade  surface  are  a contour  map  of  the  blade's  displacement  normal 
to  the  plane  of  the  hologram.  The  brightest  fringes  are  nodal  lines 
(zero  displacement)  for  the  particular  mode  of  vibration.  As  one  travels 
away  from  a nodal  line,  the  dark  fringes  represent  increasing  values  of 
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the  normal  displacement  and  are  approximately  equal  to  Integer  multiples 
of  one-quarter  wavelength  of  the  laser  light.  The  actual  value  of  the 
normal  displacement  of  the  i-th  dark  fringe  is  given  by  the  expression 
(Ref.  7): 


Wi= 


A 

2ti(Cos  01+  62) 


(3.16) 


where  A 


wavelength  of  the  laser  light  used  to  make  the  hologram 


0|  - angle  between  the  displacement  vector  of  the  plate  and  the 

line  of  sight  of  the  observer /earner a through  the  hologram 


02  - angle  between  the  displacement  vector  of  the  plate  and  the 

laser  "object"  beam  illuminating  the  plate 

- i-th  root  of  the  zero  order  Bessel  function  corresponding 
to  fringe  i. 


If  the  specimen  to  be  tested  is  positioned  such  that  its  principal  plane 
of  vibration  is  parallel  to  the  plane  of  the  hologram,  and,  in  addition, 
is  located  such  that  its  distance  from  the  hologram  is  greater  than  five 
times  its  largest  In-plane  dimension,  the  value  of  0i,  can  be  made  very 
small  and  Cos  0],  = 1.  Further,  if  the  object  beam  for  Illuminating  the 
test  specimen  is  placed  very  close  to  the  hologram  and  the  5 to  1 ratio 
(distance  of  specimen  from  hologram  to  maximum  specimen  dimension)  men- 
tioned above  is  in  effect,  then  02  can  also  be  made  very  small  and 
Cos  02  = 1.  With  these  criteria  placed  on  experimental  setup,  equation 
(3.16)  takes  the  form 


Wi  = 


A 

4 TI 
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(3.17) 


Table  3.3  gives  the  normal 
fringes  based  on  the  above 
length  for  a HeNe  laser  of 


displacement  for  the  first  twenty  interference 
simplifying  assumptions  and  a laser  wave- 
A = 6328A  = 24.914  pin  (Ref.  8). 

TABLE  3.3 


Normal  Displacement  Values  from  Holographic  Fringes  (Ref.  8) 


Fringe  Order,  i i-th  Root  of  Jq  (fij^)  Displacement  (pin) 


1 

2.405 

4.768 

2 

5.520 

10.94 

3 

8.634 

17.12 

4 

11.79 

23.37 

5 

14.93 

29.60 

6 

18.07 

35.82 

6 

21.21 

42.06 

8 

24.35 

48.29 

9 

27.49 

54.51 

10 

30.64 

60.74 

11 

33.78 

66.97 

12 

36.92 

73.20 

13 

40.06 

79.43 

14 

43.20 

85.66 

15 

46.34 

91.89 

16 

48.48 

98.11 

17 

52.62 

104.3 

18 

55.77 

110.6 

19 

58.91 

116.8 

20 

62.05 

123.0 

Utilizing  equation  (3.1),  the  experimental  normal  displacements, 
w j , at  the  locations,  Xj , j=  1,  2,  3,  ..k,  in  the  clamped-free  direc- 
tion can  be  set  equal  to  the  series  approximation  for  the  same  deforma- 
tion given  by  equation  (3.17).  In  this  case,  the  chord  location,  y, 
is  held  constant.  Doing  so,  yields 

P 

W^(Xj)  = wj  = AiXi(Xj),  j = 1.  2,  . . . k (kip)  (3,18) 

The  only  unknowns  in  equation  (3.18)  are  the  modal  weighting  coefficients,  Af. 
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Since  we  require  that  k^,  equation  (3.18)  represents  a set  of  k non- 
homogeneous  simultaneous  equations  in  p unknowns,  i.e.,  the  modal 
weighting  coefficients,  Ai,  i = 1,  2,  . . . ,p.  Equation  (3.18)  can 
be  represented  in  terms  of  matrix  notation  as 

[X]  (A)  = (W)  (3.19) 

where  [X]  - rectangular  k x p matrix  of  the  p clamped-free  beam 
mode  functions  at  the  k spanwise  locations, 

xj,  j=l,2, k. 

^[X]  has  rank  r^p) 

(A)  - column  vector  of  the  p modal  weighting  coefficients,  Ai, 
i=l,2. .p. 

(W)  - column  vector  of  the  k normal  displac'^ments , wj  , j=l,2,...k, 
obtained  at  the  locations  xj  using  holographic  inter- 
ferometry. 

For  the  case  of  k=p,  equation  (3.19)  can  be  solved  for  the  modal  weighting 
coefficients,  (A),  using  standard  matrix  inversion  techniques.  But  since 
we  are  dealing  with  data  points  determined  experimentally  and  the  result- 
ing inherent  error,  it  is  prudent  to  choose  k>p  such  that  equation  (3.19) 
represents  an  overdetermined  set  of  equations.  A "best  fit"  approxima- 
tion for  the  p weighting  coefficients,  (A)  can  then  be  found  based  on  the 
k normal  displacement  values,  (W) . This  type  of  approach  was  implemented 
using  a linear  least  square  solution  of  equation  (3.18)  as  described  in 
the  following. 
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The  k simultaneous  linear  equations  given  by  equation  (3.18)  can 
be  solved  for  the  p modal  weighting  coefficients,  (A),  for  the  case  of 
k>p  by  utilizing  the  concept  of  the  pseudoinverse  developed  by  E.  H. 

Moore  (Ref.  9)  and  R.  Penrose  (Ref.  10)  independently  of  each  other 
under  the  name  "generalized  inverse".  At  present,  the  name  "pseudo- 
inverse" is  the  more  generally  accepted  term  and  will  be  used  in  the 
following  discussion.  The  application  of  the  pseudoinverse  to  the 
solution  of  overdetermined  systems  of  linear  equations  is  dealt  with 
quite  lucidly  by  Golub  and  Reinsch  (Ref.  11). 

The  rectangular  (k  x p)  matrix,  [X],  appearing  in  equation  (3.19) 
can  be  represented  as  (Ref.  11)* 

[X]  = [U]  [X]  [V]^  (3.20) 

where  [U]  - matrix  of  the  p orthonormalized  eigenvectors  of 
T 

[X]  [X]  based  on  the  p largest  eigenvalues  of 
[X]  [X]^ 


[V]  - matrix  of  the  p orthonormalized  eigenvectors 
of  [X]’^  [X] 

[X]  - Diag  (Xj,  ....)  is  the  diagonal  matrix 

consisting  of  the  non-negative  square  roots 

T 

of  the  eigenvalues  of  [X]  [X]. 


[X]  is  arranged  such  that  X|,>  X2  > .....\p 
[U]  = [V]'^  [V]»  [V]  [V]'^  = Ip  (Ip  is  the  identity  matrix) 


* [V]  - Transpose  of  matrix  [V],  ^ 
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The  pseudoinverse  of  [X] , denoted  as  [X]  , is  obtained  from  the  ex- 
pression 

[X]"^  = [V]  [U]’^  (3.21) 


where  [U]  and  [V]  are  defined  as  before  and 


l^Y 


T-  0 0 O , 

Al 

o -J-  O 0 . 

1 

. . . o o 


[I]  [Z]  ^ = [Z]~^  [Z]  = I, 


equation  (3.21)  can  be  used  to  solve  for  the  modal  weighting 
coefficients  in  equation  (3.19). 

Pre-mult Iplylng  by  [X]^  yields: 

[X]"^  [X]  (A)  = [X]"^  (W) 

[Ip]  (A)  = [X]'^  (W) 

(A)  = [V]  [Z]“^  [U]''^  (W) 


(3.22) 

(3.23) 

(3.24) 


Equation  (3.24)  gives  the  p weighting  coefficients  based  on  the  k 
input  displacements  and  the  p clamped-free  beam  functions  evaluated 
at  the  locations  x j , j=l,2,...k.  The  nature  of  this  solution  is 
Indeed  fortuitous  for  the  following  two  reasons  (Ref.  12).  First, 
if  the  rank,  r,  of  matrix  [X]  is  r=p,  then  the  weighting  coefficient 
vector  (A)  in  equation  (3.24)  is  the  "best  fit"  in  the  sense  that  it 
is  nearest  to  the  actual  vector  (A)  in  a least  squares  sense. 

Second,  if  the  rank  of  matrix  [X]  is  r<p,  then  there  is  no  exact  solution 
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but  equation  (3.24)  still  gives  the  best  solution  in  the  sense  of 
least  squares,  i.e.,  [X] (A)  is  as  close  to  (W)  as  possible.  Appen- 
dix A gives  a simple  example  of  computing  the  pseudoinverse  of  a 3x2 
matrix. 

The  computation  of  the  psudoinverse  of  [X]  in  equation  (3.19) 
and  the  resulting  solution  for  the  modal  weighting  coefficients, 

(A),  was  carried  out  using  the  computer  program  LLSQAR  that  is  on 
permanent  file  on  the  computer  system  at  Wright-Patterson  AFB.  The 
routine  is  part  of  the  Computer  Science  Center's  International 
Mathematical  Scientific  Library  package.  A brief  description  on  the 
actual  mechanics  of  using  LLSQAR  is  given  in  Appendix  D. 

With  the  values  of  the  modal  weighting  coefficients  determined, 
equations  (3.1)  and  (3.7)  can  be  used  to  determine  the  displacement 
and  curvature,  respectively,  for  any  location,  x,  y=constant,  along  the  span 
of  the  plate. 

In  a manner  exactly  paralleling  the  preceding  development  for  plate 
deformation  based  on  a fixed-free  modal  series,  the  modal  weighting  co- 
efficients appearing  in  equation  (3.3)  for  a free-free  beam  series  can  be 
determined  by  using  the  experimental  normal  displacements,  wj , at  each  of 
k locations,  y j , j=l,2,3...k,  in  a line  across  the  chord  of  the  plate, 
in  the  free-free  direction.  In  this  case,  the  span  or  x location  of  the 
plate  will  be  held  constant.  Carrying  this  out,  yields  the  equation 

q 

Wy(y)  = Wj  = BiY^(yj)  (325) 

for  i=l,2,....q  and  j=l,2,....k  (k>q) 
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The  modal  weighting  coefficients,  Bi,  are  then  computed  using  the 
computer  algorltha  LLSQAR.  Once  this  has  been  accomplished,  the 
displacement  and  curvature  at  any  point  y,  x=constant,  can  be  computed 
via  equations  (3.3)  and  (3.9).  Section  IV  of  this  report  presents  a 
detailed  discussion  of  the  steps  Involved  in  applying  the  method  to  a 
typical  plate  vibration  mode  shape. 


3.3  Ramifications  of  Using  Simple  Beam  Eigenfunctions 

At  the  beginning  of  this  section  mention  was  made  of  the  fact 
that  while  the  clamped-free  and  free-free  beam  functions  satisfy  the 
corresponding  plate  geometric  boundary  conditions  exactly,  the  plate’s 
natural  boundary  conditions  along  its  free  edges  are  only  approximately 
satisfied.  The  ramifications  of  this  can  be  seen  by  looking  at  the 
natural  boundary  conditions  that  must  be  satisfied  along  the  tip  of  the 
cantiliver  plate  shown  in  Figure  3.1.  Satisfying  equilibrium  along  this 
‘edge  requires  that  the  shear  forces  in  the  x direction  and  the  moments 


about  the  y axis  be  zero.  Hence,  we  require  that: 

I / . 3^W\  I 


ia  L 


(3.26) 


r3^w 

\=a  = 


+ (2  - V)  — — 2 I 

3x  3y‘‘  x= 


(3.27) 


where  D Eh  „ is  the  plate  flexural  stiffness.  Now  the  clamped-free 

“ 12(l-v^) 

beam  function  given  by  equation  (3.2)  has  the  property  that  the  first 
term  in  each  of  equations  (3.26)  and  (3.27)  is  identically  equal  to  zero. 


i.e. , 


3^W  3^ 

■ 3x3  ■ »• 


(3.28) 


when  W(x)  is  approximated  by  a series  of  clamped-free  beam  functions. 
This  means  that  the  calculated  moment  and,  hence  the  strain,  e^,  at 
the  plate  free  edge  will  have  an  error  that  is  proportional  to 


6y2 


and  the  calculated  shear  will  have  an  error  that  is  proportional 


to  (2-v)i_li_^.  There  will  be  a similar  error  in  the  value  of  the  chord- 
<Sx6y 

wise  strain,  Ey,  at  the  plate's  free  edges  y=o  and  b.  Fortunately,  the 
magnitude  of  this  inherent  error  is  not  appreciable  as  will  be  seen  in 
Section  IV, 
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SECTION  IV 


INTERACTIVE  COMPUTER  PROGRAM  - HOLOCURVE 


The  technique  developed  in  Section  III  for  determining  the  displace- 
ment and  curvature  of  a vibrating  plate-like  structure  was  programmed  for 
use  on  a CDC6600  computer  via  a Tektronix  4010  interactive  computer  ter- 
minal. By  using  the  Tektronix  computer  unit,  one  has  the  advantage  of 
interacting  with  the  computer  program  while  it  is  executing  on  the 
CDC6600  computer.  In  this  way,  input  data  can  be  deleted  or  modified, 
decisions  as  to  program  direction  can  be  made,  and,  most  importantly,  the 
results  are  returned  to  the  user  immediately  in  either  tabular  or  plot 
format . 


Figure  4.1  shows  a flow  chart  of  the  computer  program  called  HOLO- 
CURVE. A listing  of  the  program  is  given  in  Appendix  B.  The  questions 
that  the  user  must  answer  while  executing  HOLOCURVE  are  listed  below  in 
order  of  appearance  on  the  Tektronix  4010  screen  along  with  possible 
answers  (in  brackets)  and  a brief  discussion  of  each  entry. 


1.  WHAT  IS  THE  MODE  NAME  (8H)?  


Input  any  letter/number  identifier  that  takes  up  eight 
spaces  or  less. 


2.  WILL  OUTPUT  BE  BASED  ON  CLAMPED-FREE  (CF)  OR  FREE-FREE  (FF) 
BEAM  SERIES?  TYPE  CF  OR  FF.  [CF  / FF] 


Program  is  henceforth  keyed  to  either  a clamped-free 
or  free-free  modal  beam  series. 
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3.  HOW  MANY  CF  [OR  FF]  BEAM  SERIES  TERMS  WILL  BE  USED?  

[Floating  point  whole  nvimber] 

This  specifies  the  value  of  p In  equation  (3.1)  for  a 
clamped-free  beam  or  of  q In  equation  (3.3)  for  a 
free-free  beam.  For  the  lower  mode  shapes  (first  and 
second  bending,  first  and  second  torsion,  lyre  or  two- 
strlpe) , a value  of  p = 5 or  Q = 5 gives  adequate  results. 

4.  INPUT  NUMBER  OF  FRINGE  V.ALUES. 

Specifies  number  of  discrete  holographic  fringe  value 
data  points  that  will  be  used. 

5.  INPUT  SPAN  [OR  CHORD]  LENGTH.  [Floating  point  number] 

For  the  clamped-free  case,  the  actual  length  of  the  span 
segment  Is  Input.  This  corresponds  to  the  beam  length,  a. 
In  equation  (3.2).  For  example,  the  distance  from  the  root 
to  the  tip  of  a compressor  blade  could  be  entered.  For  the 
free-free  case,  the  length  of  the  chord  segment  Is  Input. 
This  corresponds  to  the  beam  length,  b.  In  equations  (3.5) 
and  (3.6).  It  could,  for  example,  be  the  distance  from  the 
leading  edge  to  the  trailing  edge  of  a compressor  blade. 

6.  INPUT  CONSTANT  CHORD  (Y)  [OR  SPAN  (X)]  LOCATION.  

This  establishes  the  constant  y or  x location  on  the  plate 
or  blade  for  the  case  of  a clamped-free  or  free-free  beam 
series  approximation,  respectively.  (See  Figure  3.1) 
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7.  INPUT  PAIRS  OF  SPAN  [CHORD]  LOCATION,  FRINGE  VALUE  DATA. 

FREE  FIELD  FORMAT  WITH  EACH  VALUE  SEPARATED  BY  COMMAS. 

DATA  = [1,  -5.,  2,  -II.,  3.,  -12.,  4.1,  -8,  4.94,  0.,  6.,  9., 

7.,  16.] 

Pairs  of  span  (or  chord)  location  vs.  holographic 
fringe  value  are  input.  Seven  such  pairs  of  data  are 
input  in  the  example  shovm.  Note  that  the  fringe  values 
are  positive  or  negative  whole  numbers.  Note,  also,  there 
is  a nodal  point  at  x = 4.94  where  the  fringe  value  is 
given  as  zero.  The  HOLOCURVE  subroutine,  CONVRT,  converts 
the  fringe  values  to  normal  displacments  using  the  expression 
given  in  equation  (3.17). 

8.  WILL  DISTANCE  FROM  NEUTRAL  AXIS  VS.  SPAN  [OR  CHORD]  LOCATION 

BE  INPUT  - YES  OR  NO?  

If  "YES",  the  following  two  statements  must  be  answered: 

HOW  MANY  (SEGMENT  LOCATION /NEUTRAL  AXIS  DISTANCE)  VALUES  WILL 
BE  INPUT  (MAX.  OF  20)? [5.] 

INPUT  5.  PAIRS  OF  (SEGMENT  LOCATION,  NEUTRAL  AXIS 
DISTANCE)  VALUES. 

DATA  = [0.,  .25,  1.,  .28,  3.,  .30,  5.,  .25,  6.,  .20,  7.,  .15] 

In  order  to  compute  the  bending  strain  at  a point  on  the  plate 
or  blade,  the  distance  from  the  neutral  axis  must  be  known 
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of  a flat  plate-like  structure  having  constant  thickness, 
this  distance  is  merely  half  of  the  thickness.  In  the 
case  of  a twisted  or  curved  plate-like  body  having 
variable  thickness,  however,  the  neutral  axis  distance  is 
not  so  easily  computed  and  is  entered  in  discrete  form  as 
in  the  above  "DATA"  entry.  When  the  bending  strains  are 
computed  along  the  span  or  chord  section  being  analyzed,  a 
subroutine  (INTERP)  utilizes  a second  degree  polynomial 
interpolation  scheme  to  compute  the  neutral  axis  distance 
based  on  the  discrete  input  data.  Appendix  E.  contains 
a more  complete  description  of  subroutine  INTERP. 

9.  DO  YOU  WISH  TO  SEE  INPUT  DATA  (YES  OR  NO) ? 

A listing  of  fringe  values  and  their  location,  and,  if  used, 
neutral  axis  distance  values  and  their  location  are  displayed. 

10.  COMPUTATION  COMPLETED  FOR  [Mode  Name]  MODE,  DATA  SET  = 1.  DO 

YOU  WISH  TO  SEE  MODAL  WEIGHTING  COEFFICIENTS  FOR  THE  CF[0R  FF]  BEAM 
SERIES  - YES  OR  NO?  

An  affirmative  answer  yields  a listing  of  the  modal 
weighting  coefficients  given  in  equations  (3.1)  or  (3.3). 

11.  DISPLACEMENT  AND  CURVATURE  AS  A FUNCTION  OF  POSITION  WILL  NOW 
BE  LISTED  AND  PLOTTED  FOR  MODE  [Mode  Name]. 

HOW  MANY  DATA  POINTS  DO  YOU  WANT  (MAX.  OF  40)7  

Chord  or  span  segment  is  evenly  divided  into  number  of  data 
points  entered.  Program  then  proceeds  to  give  a listing  of 
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the  normal  displacement  as  a function  of  position  and  a 
plot  of  same.  It  then  gives  a listing  of  the  curvature 
(bending  strain)  as  a function  of  position,  and  an 
accompanying  plot. 

12.  IS  THERE  MORE  DATA  TO  BE  ANALYZED  (YES  OR  NO?)  

If  answer  la  affirmative,  cycle  of  questions  Is  repeated; 

If  negative,  HOLOCURVE  completes  Its  execution  and  user  may 
"LOG  OFF"  from  terminal.  , 


Appendix  C.  gives  a typical  example  of  the  above  statements  being 
used  during  a HOLOCURVE  execution. 
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SECTION  V 


DISPLACEMENT  AND  STRAIN  DETERMINATION  USING  HOLOCURVE 

5.1  Displacement  and  Strain  Distribution  of  a Vibrating  Plate 

In  order  to  check  the  effectiveness  and  accuracy  of  the  holographic 
data  reduction  technique  utilized  by  HOLOCURVE,  three  natural  modes  of 
vibration  of  a cantilevered  flat  plate  were  determined  using  holographic 
interferometry  and  the  resulting  fringe  data  was  input  to  HOLOCURVE  to 
determine  the  displacement  and  bending  strain  at  selected  locations  for 
each  of  the  three  vibration  modes.  The  results  from  this  analysis  were 
compared  to  those  of  an  eigenvalue  analysis  carried  out  on  a finite 
element  model  of  the  plate  using  the  multi-purpose  finite  element  com- 
puter program,  NASTRAN  (Rigid  Format  3). 

The  plate  was  made  of  6061-T6  Aluminum  and  measured  3"  by  7"  with 
a thickness  of  .19".  It  was  fixed  along  one  of  its  short  sides  by  clamp- 
ing it  between  two  massive  steel  blocks.  Time  average  holograms  were 
made  of  each  of  the  three  natural  plate  modes  of  vibration  while  the 
plate  was  being  excited  by  a 40-watt  siren.  A schematic  of  the  holo- 
graphic set-up  is  shown  in  Figure  5.1.  The  set-up  is  quite  typical  of 
the  type  used  for  doing  holographic  interferometry.  A 50mw  continuous 
wave  HeNe  laser  (Spectra-Physlcs  125A)  was  used  for  the  tests  with  the 
holograms  being  recorded  on  4"x5"  photographic  plates  (Agfa-Gevaert 
10E75). 

What  follows  are  the  results  of  using  HOLOCURVE  to  get  the  displace- 
ment and  corresponding  bending  strain  for  the  plate  vibration  modes — 
second  bending,  second  torsion,  and  the  lyre  mode. 


Schematic  of  Experimental  Set-up;  BS-beam  splitter, 
BE-beam  expander,  P-photographic  plate,  S-siren, 
beam,  0-object  beam 


5.2  Lyre  Mode  of  Vibration 


Figure  5.2  shows  the  lyre  mode  (after  the  Greek  instrument)  of  the 
cantilever  plate  obtained  using  holographic  interferometry  and  finite 
element  analysis.  Note  that  the  natural  frequencies  determined  by  each 
of  the  techniques  differ  by  only  1.2%.  In  the  figure,  the  plate  is 
clamped  along  its  bottom  edge.  The  lyre  mode,  or  a form  similar  to  it, 
is  a common  one  in  turbine  engine  compressor  blading. 

Referring  to  Figure  5.2,  let  a Cartesian  coordinate  system  originate 
at  the  lower  right-hand  corner  of  the  plate  with  the  x axis  lying  along 
the  plate's  right  edge  and  the  y axis  lying  along  the  bottom  (fixed) 
edge  of  the  plate.  Using  this  coordinate  system,  the  displacement  and 
corresponding  bending  strain  were  determined  with  respect  to  the  x coordi- 
nate (spanwlse  direction)  while  keeping  y (the  chordwlse  direction) 
constant  at  y=.5".  Note  that  the  coordinate  system  just  defined  is  the 
same  as  the  one  shown  in  Figure  3.1  of  Section  III,  Theoretical  Develop- 
ment . 

Referring  to  the  lyre  mode  shape  characterized  by  the  holographic 
fringes  in  Figure  5.2,  the  holographic  fringe  values  as  a function  of  x, 
y=.5",  were  determined  and  are  given  in  Table  5.1. 

( 

f 

TABLE  5.1  * 

Fringe  Values  vs.  x,  y=.5",  for  Lyre  Mode 
Fringe  Number,  N Span  Location,  x,  (in.) Fringe  Value,  F 


1 

.50 

1 

2 

1.25 

2 

3 

2.50 

-2 

4 

3.50 

-4 

5 

4.30 

-2 

6 

5.32 

0 

7 

6.25 

-2 

8 

6.95 

-5 
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The  fringe  values  and  span  locations  given  in  Table  5.1  are  not  the 
complete  set  of  values  occurring  on  this  spanwise  segment  of  the  plate  and, 
hence.  Table  5.1  does  not  represent  a unique  set  of  fringe  value  data  for 
the  span.  Some  general  rules  of  thumb  were  used,  however,  in  their  choice. 

1.  Fringe  value  locations  should  be  chosen  such  that  they  are 
approximately  evenly  spaced  over  the  segment  length. 

2.  The  total  number  of  fringe  values  used  must  be  greater  than  the 
number  of  beam  series  functions  input  to  HOLOCURVE.  This  will 
insure  that  full  advantage  is  taken  of  the  "smoothing"  character- 
istics of  the  linear  least  squares  subroutine,  LLSQAR. 

3.  When  possible,  specify  the  locations  of  the  nodal  lines  with 

the  exception  of  the  clamped  end  nodal  line  in  the  case  of  the 
clamped-free  beam  series  approximation.  To  do  so  in  this  latter  t 

instance  would  be  redundant  since  zero  displacement  is  one  of 

I the  geometric  boundary  conditions  of  the  clamped-free  beam  1 

I functions.  j 

4.  Finally,  regarding  the  number  of  beam  functions  specified  for  | 

HOLOCURVE,  in  the  case  of  clamped-free  beam  series  approxima- 
tion, the  number  of  beam  series  functions  used  should  be  greater 
than  or  equal  to  the  number  of  node  lines  crossing  the  span 
plus  two,  i.e. , 

Number  CF  Series  Functions  = Number  Node  Lines  + 2. 

For  the  case  of  a free-free  beam  series  approximation,  the 
number  of  beam  series  functions  used  should  be  greater  than  or 
equal  to  the  number  of  node  lines  crossing  the  chord  plus  three. 

i.e.. 

Number  FF  Series  Functions  = Number  Node  Lines  + 3. 

The  set  of  eight  fringe  values  and  their  span  locations  given  in 
Table  5.1  were  input  to  HOLOCURVE  where  they  were  used  to  compute  the  modal 
weighting  coefficients  of  a five-terra  clamped-free  beam  series.  The  com- 
puted normal  displacement  is  shown  as  the  smooth  curve  in  Figure  5.3  as  a 
function  of  the  span,  x.  Also  shown  in  the  figure  are  the  normal  displace- 
ment values  computed  from  the  finite  element  model  of  the  plate  using  NASTRAN. 

The  finite  element  values  were  normalized  to  the  maximum  displacement  at  the 
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Normal  Displacement,  W,  vs.  Span,  X,  at  Y«.5"  for  Lyre  Mode 


plate  tip  as  determined  by  HOLOCURVE.  It  is  seen  that  the  agreement 
between  the  two  methods  is  quite  good.  The  modal  weighting  coefficients 
upon  which  the  HOLOCURVE  results  are  based  are  listed  in  Table  5.2. 

TABLE  5.2 

Modal  Weighting  Coefficient  for  CF  Series 
Approximation  of  Lyre  Mode 

CF  Beam  Mode,  i Modal  Weighting  Coefficient,  AixlO^ 

8.66 
2.860 
1.886 
9.164 
1.526 

Note  that  coefficients  A^,  and  A4,  representing  the  first  and  fourth 
clamped-free  beam  bending  mode  shapes,  respectively,  are  the  dominant 
terms  in  the  series.  This  is  as  one  would  expect,  since  the  first  order 
of  deflection  of  the  mode  is  that  of  a negative  first  bending  mode  while 
the  actual  modal  displacement  is  that  of  a positive  fourth  bending  mode 
shape,  i.e. , four  nodal  lines. 

Figure  5.4  shows  the  resulting  bending  strain,  in  the  plate's 
X direction  computed  from  the  second  derivative  of  the  five-term  clamp- 
ed-free  beam  series.  As  in  the  plot  of  the  normal  displacement,  the 
individual  data  points  in  the  figure  represent  results  from  the  NASTRAN 
finite  element  analysis  of  the  plate. 


1 

2 

3 

4 

5 
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Bending  Strain, vs.  Span,  X,  at  Y=.5"  for  Lyre  Mode 


The  values  of  the  bending  strain,  computed  using  NASTRAN,  were 
normalized  to  the  HOLOCURVE  output  at  the  span  location,  x*.15". 

Again,  the  agreement  Is  seen  to  be  quite  good. 

The  displacement  and  resulting  bending  strain  were  also  determined  for 
the  lyre  mode  In  Figure  5.2  as  a function  of  the  width,  y,  for  x=2.5". 

The  holographic  fringes  upon  which  the  HOLOCURVE  computations  were 
based  are  given  In  Table  5.3. 

TABLE  5.3 

Fringe  Values  vs.  y,  x=2.5",  for  Lyre  Mode 
Fringe  Number,  N Chord  Location,  y (In.) Fringe  Value,  F 


1 

.06 

-5 

2 

.48 

-2 

3 

.80 

0 

4 

1.20 

2 

5 

1.84 

2 

6 

2.25 

0 

7 

2.53 

-2 

8 

2.96 

-5 

This  set  of  eight  fringe  values  and  chord  location  data  were  used  in 
HOLOCURVE  to  compute  the  modal  weighting  coefficients  of  a 5-term 
free-free  beam  series.  The  normal  displacement  curve  computed  from 
this  series  Is  shown  as  the  smooth  curve  in  Figure  5.5  along  with  the 
normalized  grid  point  displacements  from  the  NASTRAN  finite  element 
calculations.  Note  that  the  finite  element  computations  yield  a 
greater  value  for  the  displacement  at  the  chord  midpoint.  This  maximum 
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igure  5.6  - Bending  Strain,  ^ , vs.  Chord,  Y,  at  X=2.5"  for  Lyre  Mode 


deviation  occurred  at  mid-chord  only  because  the  finite  element  values 
were  normalized  at  the  free  edge,  y=0. 


Figure  5.6  shows  the  bending  strain,  Sy,  as  a function  of  chord 
position,  y,  for  x=2.5".  The  finite  element  data  in  Figure  5.6  were 
normalized  to  the  HOLOCURVE  strain  at  y=1.35".  As  in  the  case  of  the 
spanwise  bending  strain,  ex»  good  agreement  was  found  between  the 
HOLOCURVE  and  normalized  finite  element  results.  Table  5.4  gives  the 
modal  weighting  coefficients  from  which  the  normal  displacement  and 
bending  strain  in  the  preceding  two  figures  were  computed. 


TABLE  5.4 


Modal  Weighting  Coefficients  for  FF  Series 
Approximation  of  Lyre  Mode 


FF  Beam  Mode,  i 


Modal  Weighting  Coefficient,  BixlO 


1 

2 

3 

4 

5 


-3.45 
- .40 
-14.09 
.20 
-.04 


Note  that  coefficient,  B3,  corresponding  to  the  first  bending  mode  shape 
of  a free-free  beam  is  the  dominant  term  in  the  series  as  one  would  expect. 


5 . 3 Second  Bending  Mode  of  Vibration 

Figure  5.7  shows  the  second  bending  mode  of  vibration  for  the  canti- 
lever plate  obtained  using  holographic  interferometry  and  finite  element 
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analysis.  As  in  the  preceding  section  on  the  lyre  mode,  the  coordinate 
system  for  locating  fringe  values  and  interpreting  displacement  and 
strain  output  is  located  at  the  lower  right-hand  corner  of  the  plate. 
The  X axis  lies  along  the  plate's  right  side  and  y axis  lies  along  the 
lower  (fixed)  edge  of  the  plate. 

The  displacement  and  corresponding  bending  strain  were  determined 
with  respect  to  x (spanwise  direction)  for  y=.5"  using  HOLOCURVE  and 
holographic  fringe  values  taken  from  Figure  5.7.  The  fringe  values 
used  in  the  HOLOCURVE  program  are  given  in  Table  5.5. 


TABLE  5.5 

Fringe  Values  vs.  x,  y=.5",  for  Second  Bending  Mode 


Number,  N 

Span  Location,  x (in.) 

Fringe  V< 

1 

.44 

1 

2 

1.16 

4 

3 

2.11 

9 

4 

2.66 

11 

5 

3.88 

11 

6 

4.66 

7 

7 

5.44 

0 

8 

6.11 

-7 

9 

6.98 

-16 

This  set  of  data  was  input  to  HOLOCURVE  and  a five-term  clamped- 
free  beam  series  was  generated.  The  modal  weighting  coefficients  for 
this  beam  series  are  given  in  Table  5.6. 


TABLE  5.6 


Modal  Weighting  Coefficients  for  CF  Series 
Approximation  of  Second  Bending  Mode 


CF  Beam  Mode,  i 

Modal  Weighting  Coefficient,  AixlO^ 

1 

-.68 

2 

50.64 

3 

.83 

4 

-.24 

5 

.51 

Looking  at  Table  5.6,  it  is  evident  that  the  plate's  second  bending  mode 
shape  is  quite  similar  to  that  of  a simple  clamped-free  beam  and,  hence, 
the  modal  weighting  coefficient,  A2,  dominates  the  clamped-free  beam 
series. 


Figure  5.8  shows  the  plate's  second  bending  mode  normal  displacement 
computed  by  HOLOCURVE  as  a function  of  the  span  coordinate,  x,  for  y=.5". 
The  normal  displacement  values  determined  by  finite  element  analysis  were 
normalized  to  the  HOLOCURVE  displacement  at  the  plate's  free  edge  (x=7.0"). 
The  agreement  between  the  two  methods  of  analysis  for  this  case  are  remark- 
ably good. 

The  bending  strain.  Ex,  resulting  from  the  plate  displacement  shown 
in  Figure  5.8  is  shown  in  Figure  5.9.  As  in  the  case  of  the  lyre  mode, 
the  bending  strain  values  computed  using  the  finite  element  method  were 
normalized  to  the  HOLOCURVE  strain  value  at  x=.15".  Again,  the  agreement 
between  the  two  methods  is  seen  to  be  quite  good. 


41 


42 


It  was  not  possible  to  use  HOLOCURVE  to  compute  the  displacement 
and  bending  strain  In  the  plate's  y (chordwlse)  direction  for  the  second 
bending  mode  because  the  displacement  gradient  Is  extremely  small  In 
this  direction  being  roughly  equivalent  to  only  one  holographic  fringe. 
Hence,  the  discrete  fringe  data  necessary  for  the  HOLOCURVE  computations 
Is  nonexistent.  The  one  saving  grace  In  this  case  Is  that,  since  the 
displacement  gradient  Is  so  small  In  the  chordwlse  direction,  so,  too. 

Is  the  bending  strain  and,  hence,  the  Information  lost  would  not  be  of 
any  great  value  In  a structural  vibration  analysis. 

5.4  Second  Torsional  Mode  of  Vibration 

Figure  5.10  shows  the  second  torsional  mode  of  vibration  for  the 
cantilever  plate  obtained  using  both  holographic  Interferometry  and  the 
finite  element  method.  This  mode  of  vibration  was  chosen  as  a HOLOCURVE 
application  because  It  points  out  a characteristic  of  torsional  vibration 
that  Impacts  HOLOCURVE  results  for  chordwlse  calculations.  This  point 
will  be  addressed  In  detail  below. 

Using  the  same  coordinate  system  as  In  the  preceding  two  vibration 
modes,  the  mode  displacement  and  bending  strain  were  computed  by  HOLO- 
CURVE as  a function  of  the  span  position,  x,  for  y=.5".  The  holographic 
fringes  used  for  the  computations  were  taken  from  the  inter ferogram  in 
Figure  5.10  and  are  shown  in  Table  5.7. 
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TABLE  5.7 

Fringe  Values  vs.  x,  y=.5",  for  Second  Torsion  Mode 
Fringe  Number,  N Span  Location,  x (In.) Fringe  Value,  F 


1 

1.00 

-5 

2 

2.00 

-11 

3 

3.00 

-12 

4 

4.00 

-8 

5 

4.94 

0 

6 

6.00 

9 

7 

7.00 

16 

Entering  this  set  of  data  Into  the  HOLOCURVE,  the  modal  weighting 
coefficients  were  computed  for  a five-term  clamped-free  beam  series  and 
are  shown  in  Table  5.8. 


TABLE  5.8 

Modal  Weighting  Coefficients  for  CF  Series 
Approximation  of  Second  Torsion  Mode 


CF  Beam  Mode,  1 


Modal  Weighting  Coefficient,  A-jxlO” 


1 

2 

3 

4 

5 


11.72 

-50.84 

-12.54 

-.46 

-1.49 


This  five-term  series  was  used  to  compute  the  normal  displacement,  W, 


J 
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and  corresponding  bending  strain,  e^,  which  are  shown  as  the  smooth 
curves  in  Figures  5.11  and  5.12,  respectively.  The  HOLOCURVE  and 
finite  element  results  are  again  in  good  agreement.  As  in  the  pre- 
ceding two  cases,  the  finite  element  data  was  normalized  to  the  HOLO- 
CURVE displacement  with  respect  to  the  plate's  free  end  and  to  the 


HOLOCURVE  bending  strain  at  x=.15". 


An  analysis  was  also  carried  out  using  HOLOCURVE  to  determine  the 
displacement  and  bending  strain  in  the  plate's  y (chordwise)  direction 
with  x=1.0".  Table  5.9  lists  the  holographic  fringe  data  used  in  the 
analysis. 

TABLE  5.9 

Fringe  Values  vs.  y,  x=1.0",  for  Second  Torsion  Mode 
Fringe  Number,  N Chord  Location,  Y (in.) Fringe  Value,  F 


1 

.08 

-7 

2 

.30 

-6 

3 

.50 

-5 

4 

.70 

-4 

5 

.90 

-3 

6 

1.10 

-2 

7 

1.32 

-1 

8 

1.47 

0 

9 

1.62 

1 

10 

1.83 

2 

11 

2.02 

3 

12 

2.22 

4 

13 

2.43 

4 

14 

2.62 

6 

15 

2.83 

7 

16 

3.00 

8 

This  set  of  fringe  data  differs  from  the  others  used  up  to  this  point  in 
that  all  of  the  fringes  and  their  corresponding  chord  locations  at  x=1.0" 
were  used. 


47 


Bending  Strain,  fc^,  vs.  Span,  X,  at  Y-.5"  for  Second  Torsion  Mode 


An  eight-term  free-free  beam  series  based  on  the  above  data  was 
generated  by  HOLOCURVE.  The  resulting  modal  weighting  coefficients 


are  listed  in  Table  5.10. 


TABLE  5.10 

Modal  Weighting  Coefficients  for  FF  Series 
Approximation  of  Second  Torsion  Mode 


FF  Beam  Mode,  i 


Modal  Weighting  Coefficient,  BjxlO 


-26.63 


Looking  at  Table  5.10,  it  is  seen  that  the  modal  weighting  coefficient, 
B2,  representing  rigid  body  rotation,  dominates  the  free-free  beam  series. 
This  fact  is  further  underlined  by  looking  at  the  normal  displacement  as 
a function  of  y,  x=1.0",  computed  from  using  the  eight-term  beam  series 
shown  in  Figure  5.13.  The  normal  displacement  is  a sloping  straight  line 
passing  through  zero  at  mid-chord  (y=1.5").  Further,  the  NASTRAN  generated 
displacement  follows  the  HOLOCURVE  result  quite  closely. 


The  question  now  arises  that  if  the  normal  displacement  is  a constant 
or,  at  most,  a linear  function  of  the  spatial  variable,  y,  (see  equations 
3. A and  3.5),  what  meaning  does  it  have  to  take  the  second  derivative  of 
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the  free-free  beam  series  in  order  to  obtain  the  bending  strain?  The 
answer,  of  course,  is  that  the  bending  strain  output  in  this  case  is 
meaningless.  This  is  because  the  strain  is  computed  from  the  insigni- 
ficant terms  in  the  free-free  beam  series.  This  is  borne  out  by  a look 
at  the  bending  strain,  Cy,  computed  by  HOLOCURVE  and  based  on  the  modal 
weighting  coefficients  in  Table  5.10  as  shown  in  Figure  5.14.  The  HOLO- 
CURVE strain  values  oscillate  about  zero  as  one  travels  in  the  y direc- 
tion and  do  not  agree  with  the  NASTRAN  generated  strains.  Physically 
speaking,  this  anomalous  result  is  due  to  tne  fact  that  it  is  the  shear- 
ing stress  or  strain  that  is  a maximum  in  the  chordwlse  direction  for  a 
purely  torsional  mode  shape,  thus,  making  the  principal  stresses  (and 
curvatures)  a minimum.  NASTRAN  gives  reasonable,  albeit  small,  strain 
values  because  it  relies  on  a two-dimensional  plane  stress  theory  while 
HOLOCURVE  utilizes  a simple  beam  theory.  The  saving  grace  in  this 
instance  is  that  the  strain  is  small.  In  other  words,  along  the  plate 
directions  where  the  bending  strain  amplitude  is  moderate-to-large, 
HOLOCURVE  is  readily  applicable. 

As  a precaution  against  use  of  strain  values  based  on  Insignificant 
modal  weighting  coefficients,  HOLOCURVE  determines  the  absolute  value  of 
all  of  the  modal  weighting  coefficients  in  the  beam  series.  If  the 
absolute  value  of  the  largest  coefficient  is  ten  times  greater  than  that 
of  the  next  coefficient,  a message  is  transmitted  to  the  program  operator. 
Further,  if  the  largest  coefficient  corresponds  to  the  translational  or 
rotational  series  component  in  the  free-free  beam  series,  a message  to 
this  effort  is  transmitted  to  the  user. 
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DISCUSSION  AND  CONCLUSION 


The  work  described  herein  has  presented  an  analytical  method  for 
determining  the  strain  distribution  over  the  surface  of  a plate-like 
structure  based  on  the  interference  fringes  resulting  from  time  average 
holography.  The  technique  has  the  advantage  that  it  is  easily  programmed 
on  a digital  computer  and  that  the  holographic  fringe  data  used  in  the 
solution  does  not  require  any  special  treatment — neither  electronic  nor 
optical. 


Within  the  governing  bounds  of  classical  plate  theory,  the  method 
was  shown  to  agree  well  with  finite  element  analysis  as  evidenced  by  the 
examples  in  Section  IV.  Caution  should  be  used,  however,  when  applying 
the  method  to  plate  deformations  that  exhibit  a high  degree  of  transla- 
tional or  rotational  motion  as  the  resulting  curvature  and  corresponding 
strain  may  be  circumspect.  The  computer  program,  HOLOCURVE,  has  a warning 
to  this  effect  if  one  of  the  translational  or  rotational  modal  weighting 
coefficients  is  ten  times  greater  than  the  next  largest  coefficient  in  the 
series  approximation. 


Future  work  in  this  vein  of  the  current  effort  should  address  the 
effect  of  curved  surfaces  on  the  accuracy  of  the  computed  strain  values. 
Also,  plate-like  structures  when  the  neutral  axis  and  cross-sectional 
center  of  gravity  do  not  coincide  should  be  studied. 

In  conclusion,  it  is  felt  that  the  technique  embodied  by  the  computer 
program,  HOLOCURVE,  will  provide  a useful  tool  for  analyzing  the  surface 
strain  distribution  of  vibration  mode  shapes  of  plate-like  structures  such 
as  turbine  engine  blading. 

54 


Given  a k x p matrix.  A,  one  can  compute  its  psudoinverse  as 


follows: 


1.  Given: 


1 1 


2 0 


0 2 


k = 3,  p = 2 


(A-1) 


2,  We  can  then  compute  A^A. 


’l  2 o' 
10  2 


0 3 


(A-2) 


3.  Compute  the  p non-negative  square  roots  of  the  eigenvalues, 
T 

s,  of  A A to  form  the  diagonal  matrix,  Z. 


5-s  1 


s - 10s  + 24  = 


(A-3) 


Si  = 6,  s =4 


4^  0 


0 2 


(A-4) 


R 

ti 


ll 


4.  Compute  the  orthonormalized  eigenvectors,  of  A A to  form 

the  matrix  Vjq. 


(A^A  - Si  I)  Vi  = 0,  i = 1,  2 


(A-5) 


(A-6) 


(A-7) 


(A-8) 


Similarly,  for  S2  = 4: 


Hence,  from  equations  (A-8)  and  (A-9) , we  have: 


(A-9) 


(A-10) 
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5.  Compute  the  p orthonormallzed  eigenvectors,  U^jj,  based  on  the 

T 

p largest  eigenvalues  of  AA  to  form  the  matrix  U. 


T 

The  eigenvalues  of  AA  can  be  found  to  be: 

S]^  = 6 

52  = 4 

53  = 0 

The  two  orthonormallzed  eigenvectors  based  on  the  largest  two 
eigenvalues  given  In  equation  (A-12)  can  now  be  determined. 


For  SI  = 6, 


(A-11) 


(A-12) 


(A-13) 


(A-14) 


(A-15) 


7.  The  pseudoinverse  of  A is  given  as: 


A+  = Vn  Z ^ UjJ 


Hence : 


1 1 


1 -1 


1 0 


0 1 
2 


(A-17) 


(A-18) 


It  can  be  easily  shown  that  equation  (A-18)  will 
satisfy  the  conditions  for  the  pseudoinverse  one  of 
which  is: 


(A-19) 
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4.5 


(MTA  SET 


DO  YOU  UISH  TO  SEE  INPUT  DATA  <VES  OR  NO)?. 


r 


I.  ■y»w»wpii> 
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(0 

'H 


%9 
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(L 
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IO*4M<PCDOSU> 

I 

H A 11  H M M n 
U.U.IiLli.li.U.U. 


•^CMfO^^SDN 

nnns- 
^ ^ ^ ^ 
xxxxxxx 


FOR  MODE  2 FLEX 


OJRUATURE  FOR  NODE  2 FLEX 


APPENDIX  D 
LLSQAR  DESCRIPTION 


AD-A046  153 


UNCLASSIFIED 


AIR  FORCE  AERO  PROPULSION  LAB  WRI6HT-PATTERS0N  AFB  OHIO  F/6  20/11 
QUANTITATIVE  DISPLACEMENT  AND  STRAIN  DISTRIBUTION  OF  VIBRATING  — ETC(U) 
JUL  77  J C HACBAIN 


1 


AFAPL-TR-77-44 


NL 


2cf2 

A0A048IS3 


END 

DATE 


12-77 
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SUBROUTINF  LL<:0«R  ( ft  «B»MtN4«NR  1 1 4 « IBt  lOG)  (VKARFA  « 1FR> 


I ! 
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c 
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llsoar 

FUNCTION 

USA6F 

PftPAMETEBS 


M 

NA 

NR 

lA 

IB 

lORT 

mkarea 

lEP 


PRECISION 

REQO,  JNSI.  ROUTINFS 
LANGUAGE 


LLSOOOlO 

LLSUOU^O 

•library  I— 

LLSQ0040 

- LEAST  SQUARES  SOLUTION  OF  .OVERUE lERMINEO  LLSQ0050 

SYSTEM  OF  LINEAR  EQUATIONS  LLS00060 

• CALL  LLSOAR(A«P*M*NA«NdtlAtlb*It}GTtkiKAHEA*IER)LLSu0070 

- THE  COEFFICIENT  MATRIX  OF  THF  EUUAIION  LLSOOU80 

•Xbh*  MHERE  a is  M X NA  WITH  M GREATER.  THAN  LLSQ0090 
OR  EOUAL  to  na.  INHUI  a is  RtPLACtO  HY  LLSUulOO 

THE  PSEUDO-INVERSE  OF  A LLSUOllO 

- MATRIX  OF  THE  RIGHT  HAND  SIOE  OF  THE  EUUAIION  LLSUHliiO 

AXsBt  WHERE  R IS  H X NR.  THF  NA  X NH  LLSUU130 

SOLUTION  X OVERWRITES  «.  LLSOOUO 

- NUMBER  OF  ROWS  IN  A AND  B.  LLSOOlSO 

- NUmhER  of  columns  in  matrix  a LLSUUlbO 

- FlUMHEK  OF  COLUMNS  IN  HaTRIx  H LLS00170 

- MOW  DIMENSION  OF  A IN  THE  CALLING  LLSUOISO 

RR06HAM.  LLSU0I90 

- ROW  DIMENSION  OF  H In  THE  CALLING  LLSQO^OO 

PROGRAM.  LLSU0210 

- The  ELEMENTS  OF  A ARE  ASSUMED  TO  RE  CORRECT  LLSU022b 

TO  I06T  SIGNIFICANT  DIGITS.  IDGT  IS  AN  LLSU0230 

INPUT  Parameter.  LLSOU2AO 

- WORK  area  of  OIMEMSIOn  greater  than  OR  EOUAL  LLSu0250 

TO  NAINa^A).  LLSa02b0 

- ERROR  PARamETFR  LLSU0270 

terminal  error  = 12H  ♦ N LLSU02BO 

N r I INDICATES  That  INPUT  A IS  A 2ERO  LLSO0290 
MATRIX.  LLSUOiOO 

- single  LLSUC310 

- LPSO0RtLSVALR»UEPTST. VSORTm  LLSU0320 

- fortran  LLSO>j330 

LLSU03A0 


CALL  LLSQAR  (A.B,M,NA,NB,IA,IB,IDGT,UKAR£A,IER) 

Purpose 

LLSQAR  performs  Che  least  squares  solution  of  an  overdetermined  system  of  linear  equations. 
Given  the  M by  NA  matrix  A and  Che  M by  NB  matrix  B,  where  M la  not  less  than  NA.  this  routine 
finds  Che  NA  by  NB  matrix  X such  that  ] |x| and  | |aX-B| - minimum,  whara  | | * 1 i*  iBc 

Euclidean  norm  of  a matrix.  That  la. 


AX-B 


r “ 

NB 

/ ^ 

■ 

z 

Li-i 

j-1 

Vk-1 

^lk*kj 


-.)■] 


1/2 


Algorithm 

The  routine  first  calculates  the  pseudo- Inverse  of  the  matrix  A calling  IMSL  routine  LPSDOR. 

A^  A v(S^)U^.  (See  the  documentation  of  LPSDOR.)  V la  an  NA  by  NA  orthogonal  matrix,  S^ 

T 

la  an  NA  by  NA  diagonal  matrix  and  U la  an  NA  by  M matrix  with  orthogonal  rows.  Tha  solution 

^ T 

of  the  equations  la  computed  by  V(S  )U  B. 
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PronraiMiilnK  Notes 

1.  In  the  calling  program,  both  A and  B are  two-dimensional  arrays.  The  value  of  the  first 
dimension  of  A Is  IA(1A  not  less  than  M) , and  the  value  of  the  first  dimension  of  B Is 
IB (IB  not  less  than  M). 

2.  The  work  area,  WKAREA,  must  be  a one-dlmenslonal  array  of  length  at  least  NA(NA-t-4). 

3.  The  input  matrix  A Is  replaced  by  the  pseudo- Inverse  of  A. 

4.  The  NA  by  NB  solution  X overwrites  Input  matrix  B. 

Accuracy 

The  elements  of  A are  assumed  to  be  correct  to  IDGT  significant  digits.  The  routine  computes 

the  pseudo-inverse  of  A = USV  ; A differs  from  A In  the  Euclidean  norm  by  less  than  d*i?1^NA  * 
~IDGT 

(max  a.,)  * 10  , did  of  all  matrices  which  differ  from  A by  less  than  d in  the  Euclidean 

l.J  _ 

norm,  A has  tolnimal  rank. 

DIMENSION  A(6,3),B(7,4),WKAREA(60) 

Input: 


33.0 

16.0 

72.0 

-24.0 

-10.0 

-57.0 

- 8.0 

- 4.0 

-17.0 

33.0 

16.0 

72.0 

-24.0 

-10.0 

-57.0 

- 8.0 

- 4.0 

-17.0 

1.0 

0.0 

0.0 

-359.0 

0.0 

1.0 

0.0 

281.0 

0.0 

0.0 

1.0 

85.0 

1.0 

0.0 

0.0 

-359.0 

0.0 

1.0 

0.0 

281.0 

0.0 

0.0 

1.0 

85.0 

X 

X 

X 

X 

CALL  LLSQAR  (A, B,M, NA, NB, lA, IB. IDGT, WKAREA, lER) 
Output: 


“-9. 

- 8. 
_ 2. 


-9.6667  -2.6667  -32,000 

8.0000  2.5000  25.500 

, 2.6667  .66667  9.0000 


1.0000"] 

-2.0000 

-5.OOO0J 


Note:  X denotes  an  element  which  Is  not  used  by  LLSQAR. 


SUBROUTINE  LPSOOR  I AtMtNt lAtAiNVt IDGTiWKAREA* lER) 


r 

C-LPSDOR S 


r 

c 

c 

c 

c 

c 

c 

c 

c 

c 

r 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


FUNCTION 

usage 

PARAMETERS  A 
M 
N 
lA 


AINV 

lOGT 


MKAREA 

lEP 


PRECISION 

REOO.  IMSL  ROUTINES 
LANGUAGE 


library  

- PSEUOO-INVERSE  OF  A MATRIX. 

- CALL  LPSOOR  (AtMtNflAtAlNVtlOGTtMKAREAtlER) 

- A GIVEN  matrix  of  SIZE  M X N 
• NUmRER  of  ROMS  IN  MATRIX  A 

- number  of  columns  IN  Matrix  a 

- ROW  OIHENSION  OF  A AND  AINV  IN  CALLING 

PROGRAM  (lA  MUST  BE  GREATER  THAN  OR  EQUAL 
TO  M), 

- THE  TRANSPOSE  OF  THE  PSEUDO-INVERSE  OF  A 

IS  STORED  IN  THE  M X N MATRIX  AINV. 

- The  elements  of  a are  assumed  to  he  correct 

TO  I06T  DECIMAL  PLACES.  lOGT  IS  AN  INPUT 

parameter. 

- WORK  area  of  DIMENSION  AT  LEAST  EQUAL  TO 

NIN«4)  . 

- ERROR  parameter 
TERMINAL  ERROR  > 1?8  « N 

N = I INDICATES  THAT  A IS  THE  ZERO  MATRIX 
TO  lOGT  ACCURACY! IN  NORM). 

- single 

- vsortm*lsvalr*uertst 

- FORTRAN 


LPSOOOlO 

LPSU0020 

LPS00030 

LPSOOOAO 

LPS00050 

LPS00060 

LPS00070 

LPS00060 

LPS00090 

LPSOOIOO 

LPSOUllO 

LPSooiao 

LPSi)0130 

LPS001*0 

LPSDOIS.' 

LPSD0160 

LPSO0I70 

LPSUOlHO 

LPSooisr 

LPSUU200 

LPSD0210 

LPSD0220 

LPS002J0 

LPSD024G 

LPS00250 

LPSD0260 

LPSD0270 


CALL  LFSDOR  (A,M,N,IA,AINV,IDGT,WKAR£A,1ER) 

Purpose 

LPSDOR  calculates  the  pseudo- Inverse  of  the  M by  N matrix  A with  M not  less  than  N. 

Algorithm 

The  routine  calculates  the  singular  value  decomposition  of  A by  calling  IMSL  routine  LSVAIH. 

T 

A-USV  , where  U Is  an  M by  N matrix  with  orthononoallzed  columns,  S is  an  N by  N diagonal 
matrix  containing  the  singular  values  of  A,  and  V Is  an  N by  N orthogonal  matrix.  The  pseudo- 
inverse Is  given  by  V(S  ^)U^. 

Programming  Notes 

1.  In  the  main  program,  both  A and  AINV  are  assumed  to  be  of  dimension  M by  N (lA  not  less 
than  M).  Thus  the  transpose  of  the  pseudo-inverse  of  A Is  stored  In  AINV.  AINT  may  occupy 
the  same  storage  as  does  A,  or  completely  separate  storage. 

2.  The  work  area,  HKAREA,  must  be  a vector  of  length  at  least  N(Nl-A). 

3.  If  A Is  the  zero  matrix,  a terminal  error  message  Is  printed,  and  tho  zero  pscudu- Inverse 
is  not  stored  In  AINV.  Instead  It  will  contain  U. 

Accuracy 

The  elements  of  A arc  assumed  correct  to  IDGT  significant  digits.  I'he  routine  computes  tho 

exact  pseudo-inverse,  without  any  roundoff  error,  to  a matrix  A which  differs  by  A In  the 

••IDGT 

Euclidean  norm  by  less  than  10  , and  which  has  minimal  rank. 


3 

J 


I 


S 
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APPENDIX  E 


ij 


INTERP  DESCRIPTION 


Identification; 


INTKRP  Aitken's  Kth  Degree  Polynomial  Interpolation  of  Tabular  Data 
CDC  6600  FORTRAN  Subroutine  Subprogram 

Purpose: 

Given  an  ascending  or  descending  table  X(I)  of  independent  variables 
and  a corresponding  table  'i(l)  of  dependent  variables,  subroutine  INTthl' 
computes  a Kth  degree  polynomial  interpolated  value  iO  for  a given  abscissa 
XO.  If  the  argtunent  falls  outside  the  range  of  tabular  values,  extrapola- 
tion is  performed. 

Control: 

DIMENSION  X(N),  r(N),  TEMP(2{K+l)) 

CALL  INTERP  (X,Y,N,K,XO,YC ,TEI>IP,IER) 

where ; 

X(l)  List  of  values  of  the  Independent  variable  in  either 
ascending  or  descending  order 

Y(I)  List  of  the  corresponding  values  of  the  dependent  variable 
N Number  of  X,  Y pairs 

K Degree  of  the  interpolating  polynomial  (K  less  than  N) 

XO  Point  at  which  interpolation  is  requested 
YO  The  computed  interpolated  value 

TEMP  A one-dimensional  array  of  2*(K+l)  words  of  temporary  storage 

lER  = 0 Interpolation  successfully  performed 

= 1 Extrapolation  successfully  performed 

= 2 Either  degree  of  interpolation  out  of  range  or  two  identical 
independent  variables . 

Method: 

Repeated  bisection  of  the  index  (subscript)  of  the  list  of  values  of  X 
is  carried  out  until  the  argument  is  isolated  between  two  consecutive  values. 
Then  Aitken's  recursion  of  linear  interpolation  is  performed  using  an  equal 
number  of  points  on  either  side  of  XO  and,  for  even  degrees,  one  additional 
point  nearest  XO. 

Remarks ; 

The  user  is  cautioned  that  increasing  the  degree  of  the  fitted  poly- 
nomial does  not  necessarily  increase  the  "accuracy"  of  the  interpolation. 
Indeed,  this  may  introduce  spurious  oscillations  into  the  interpolatovi  values 
since  the  higher  the  degree  of  a polynomial,  the  more  relative  maxima  and 
minima  it  can  have.  In  particular,  if  the  tabular  data  involves  appivciabio 
"noise",  then  use  of  a technique  wnich  "smooths"  rather  than  one  which  "fits" 
(as  does  INTERP)  may  be  more  appropriate. 

The  user  is  also  ••neouraced  to  test  to  sor  if  IKI-;  = if  taorc  is  :u.y 
.iU  ’jti  or  .•.•nccrr.i".,’  the  .:ii  ;u. -lii.s:;  of  ttie  indcpundout  vai'i-ibleu . 

Storage ; 

INTERP  uses  23^g  words  of  storage  and  no  common. 

References; 

F.  B.  Hildebrand,  Introduction  to  Numerical  Analysis,  McGraw-Hill, 

New  York/Toronto/London,  19'.>ti,  pp.  U9-'i0. 
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